Lattice QCD at the physical point: light quark masses 



S. Durr a > c , Z. Fodor a < c ' d , C. Hoelbling a ' b , S.D. Katz d - a , S. Krieg a ' c , T. Kurth a , 
L. Lellouch b , T. Lippert c > a , K.K. Szabo a , G. Vulvert b 

a Bergische Universitdt Wuppertal, Gaussstr. 20, D-4-2119 Wuppertal, Germany. 
b Centre de Physique Theorique 1 , Case 907, Campus de Luminy, F-13288 Marseille, France. 
c Jiilich Supercomputing Centre, Forschungszentrum Jiilich, D-52425 Julich, Germany. 
d Institute for Theoretical Physics, Eotvos University, H-1117 Budapest, Hungary. 

Budapest-Marseille- Wuppertal Collaboration 



Abstract 

Ordinary matter is described by six fundamental parameters: three couplings 
(gravitational, electromagnetic and strong) and three masses: the electron's 
(m e ) and those of the up (m u ) and down (rrid) quarks. An additional mass 
enters through quantum fluctuations: the strange quark mass (m s ). The three 
couplings and m e are known with an accuracy of better than a few per mil. 
Despite their importance, m u , (their average m u( t) and m s are far less accu- 
rately known. Here we determine them with a precision below 2% by perform- 
ing ab initio lattice quantum chromodynamics (QCD) calculations, in which all 
systematics are controlled. We use pion and quark masses down to (and even 
below) their physical values, lattice sizes of up to 6 fm, and five lattice spac- 
ings to extrapolate to continuum spacetime. All necessary renormalizations are 
performed nonperturbatively. 



The masses of the up, down and strange quarks cannot be measured using 
standard experimental methods. The strong interaction confines quarks within 
hadrons (e.g. protons) in such a way that a single quark cannot be isolated. 
Moreover, the strength of the interaction is such that the mass of a hadron 
is not the simple sum of the masses of the quarks it contains. Rather it is 
provided by complicated nonperturbative dynamics (e.g. [H). This confinement 
mechanism is the low energy counterpart of the strong interaction's asymptotic 
freedom @, 01 ; by which the interactions between quarks and gluons weaken as 
their relative momenta are increased. 

Interestingly enough, the experimental data for m u , and m s has been 
available for about sixty years (the pion and kaon were discovered in the late 
1940's and the proton already 30 years before). Even the theory of the strong 
interaction, QCD, which-in principle-completely describes bound states of light 
quarks, has been known for almost four decades The fact that such a fun- 
damental question has remained poorly answered despite the available exper- 
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imental and theoretical knowledge is related to the computational difficulties 
one encounters when trying to solve the underlying equations in the domain of 
interest. The only known systematic technique to solve them is lattice QCD 
0, Several decades of theoretical, algorithmic and hardware development 
have been necessary to reach the level at which the light quark masses can be 
determined reliably. This determination is the goal of the present paper. 

For many years calculations were done in the quenched approximation. Al- 
though this approach omits the most computationally demanding part of a full 
QCD calculation -the quark determinant obtained after integrating over the 
fermion fields- a controlled determination of the strange quark mass in this ap- 
proximation (with m u =md=m s equal to about half the physical m s ) took about 
20 years Q- Moreover, the physics of the u and d quarks remained inaccessi- 
ble, because the quenched approximation, an uncontrolled truncation of QCD, 
distorts the small quark mass behavior [H, 0| . 

A very important step in the determination of light quark masses was made 
with the inclusion of u and d sea quark effects (N f =2) [lfj-llj. But even there, 
physical m u d remained elusive, this time for algorithmic reasons. A first break- 
through was made by the MILC collaboration [TBI, which used an Nf=2+1 
staggered fermion formulation to include strange sea quark effects, pushing cal- 
culations to smaller light quark masses, finer lattices and larger volumes. Up- 
dates from calculations with sea pion masses down to 258 MeV (given by the 
RMS average over taste partners for a=0.06 fm-their lightest valence pion is 
177 MeV at a=0.09 fm) and on even finer lattices are presented in (m Il7j]. 
On a subset of the MILC configurations, the HPQCD collaboration has ob- 
tained indirectly m s and m u d via the m c /m s ratio IH 1^ . Due to their use 



of quenched and partially quenched charmed and strange quarks with a non- 
unitary staggered formalism and to their error estimates on the input quantities 
that they use (m c and ri), this work does not fulfill our conditions for a con- 
trolled ab initio calculation (see below). Recently, also ETMC (Nf=2) [2(| and 
RBC-UKQCD (N f =2+1) HU have presented results with M w >270 MeV and 
significantly larger error bars. All Nf =2+1 results for m s and m u d (except 
for those of the very recent 21]) were combined into world averages in [22[, 
which also reviews Nf=2 and non-lattice results. Our results are in complete 
agreement with these averages, albeit with uncertainties smaller by more than 
a factor of 5. 

A second breakthrough came recently when it was shown that improvements 
in algorithms 23, 13] allowed the use of Wilson and domain wall fermions, which 
are free from the open questions of the rooted staggered approach, for ab initio 
calculations (e.g. [l|, |25j) and even for reaching physically light m u d, albeit in 
small volumes and at a single lattice spacing [26[. 

All previous lattice calculations of m u d and m s have neglected one or more of 
the six ingredients which we believe are most important for a full and controlled 
calculation: 

1. The inclusion of the up (u), down (d) and strange (s) quarks in the 
fermion determinant with an exact algorithm (i.e. with no integration error) 
and with an action whose universality class is QCD. Rooted staggered fermions 
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provide a numerically efficient way to investigate nonperturbative QCD. How- 
ever, this discretization is neither local nor unitary for a>0, makin g it difficult 
to show that it leads to QCD in the continuum limit (please see [16[ for an- 
other point of view). Here we use, instead, Nf =2+1 Wilson fermions with local 
improvement terms which do not affect the continuum limit. 

2. Controlled interpolations and extrapolations of the results to physical 
quark masses. Practically it means reaching pion masses as small as 200 MeV 
(clearly the value depends on the problem and on the required accuracy) or 
most preferably simulating at the physical mass point itself. At three of our 
lattice spacings we use physical (or even smaller) light quark masses. 

3. Large volumes to guarantee small finite-size effects. Our finite volume 
corrections are tiny (we use volumes up to 6 fm). Nevertheless we correct for 
them. 

4- Controlled extrapolations to the continuum limit. This requires that calcu- 
lations be performed at no less than three values of the lattice spacing, to check 
whether the scaling region is reached. We use five lattice spacings between 0.116 
and 0.054 fm, thereby gaining full control on the continuum extrapolation. 

5. Nonperturbative treatment in all steps. We obtain our primary results 
{m u d and m s in the RI scheme at 4 GeV) in a completely nonperturbative 
manner. In particular, we eliminate all truncation errors associated with the 
often used perturbative renormalization. 

6. Input parameters. The parameters of the theory (scale and quark masses) 
should be fixed with well measured observables whose error bars are undisputed 
and whose connection to experiment is transparent and contains no hidden 
assumptions. To that end we use M n , Mk and Mq exclusively. The influence 
of their error bars is negligible on our final uncertainties. Taking instead derived 
quantities, like m c and n as is done in 18, while fine in principle, can be 



problematic in practice. The error assigned to the input quantity m c in |19j is 
smaller by a factor 13 than that of the necessarily conservative Particle Data 
Group value j27|. Similarly, due to the difficulties in estimating its systematic 
uncertainty, the continuum value of r\ (and the related r ) obtained by different 
groups shows significant differences (e.g. 2.3 combined sigma between [l6[ and 



2111)- 

In this paper we determine m uc i and m s , while fulfilling all of the above 
conditions. This determination requires two, apparently straightforward, calcu- 
lations. First we compute hadron masses for tuning the quark masses to their 
physical values. Then we determine the renormalization constant to convert the 
bare quark masses to finite quantities in the continuum limit. 

We now list the most important steps of our work: 

(i) Production of the Nf=2+1 gauge field ensembles. We use a Symanzik 
improved gauge action and 2-level HEX [28j-|30j smeared clover fermions, with 
m s held close to its physical value. Gauge field configurations for 47 different 
values of the parameters ((3=6/ g 2 , am u d and am s ) were produced (c.f. Fig. Q] 
for our M,r<400MeV Nf = 2 + 1 data). 

We used five lattice spacings (owO.116, 0.093, 0.077, 0.065 and 0.054 fm), 
which are the basis for the continuum extrapolation. As we will see, the dif- 
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Figure 1: Summary of our simulation points. The pion masses and the spatial sizes of the 
lattices are shown for our five lattice spacings. The percentage labels indicate regions, in which 
the expected finite volume effect [3f| on M T is larger than 1%, 0.3% and 0.1%, respectively. 
This effect is smaller than about 0.5% for all of our runs and, as described, we corrected for 
it. Error bars arc statistical. 



ference between the results obtained on the finest lattice and those in the con- 
tinuum limit is ~3%, whereas between those of the coarsest lattice and the 
continuum limit is ~10%. 

At two pion mass points we carried out detailed finite V analyses, which 
give us a full understanding of the finite V corrections, as well as their M„ 
dependence [32j . In our calculations, we have M n L>A and/or L>5fm, so that 
the limit V— >oo can be taken safely. The difference between the results obtained 
directly on our large lattices and those in the V— >-oo limit is below the five per 
mil level. Furthermore, for M„.<200 MeV, which is most relevant for our final 
result, these corrections are even smaller, namely on the one per mil level (see 

Fig.©. 

In our calculations M„ ranges from «380 down to «120MeV (for three of 
the five lattice spacings we tuned M T to the vicinity of 135 MeV and for the 
two finest lattices, the smallest M v are around 180 and 220 MeV, respectively). 
Bracketing the physical mass point allows us to circumvent potentially trou- 
blesome chiral extrapolations. We perform calculations with m s values slightly 
below and above the physical mass, allowing a straightforward interpolation. 

(ii) Hadron and bare quark mass calculations. The pion and kaon masses 
are used to fix m u d and m 8 respectively, with Mq providing the overall scale. 
We take M,r~135 MeV, M K ~A95 MeV and Afr2~1672 MeV as input parameters 
[32l | . The calculation of hadron masses and the "mass independent scale setting" 
follows that of Tj. All three hadron masses receive finite volume corrections, 
falling off exponentially with M n L [3^]. Even though these corrections are tiny, 
they are corrected for [32j]. In addition to the hadron masses, the unrenormalized 
partially conserved axial current (PCAC) quark masses are determined. 
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(Hi) Renormalization of the bare quark masses. In addition to the PCAC 
masses discussed above, the bare m u d and m s in the Lagrangian also provide a 
measure of the quark masses used in our simulations. Once suitably renormal- 
ized, these two definitions yield quark masses which must agree in the continuum 
limit. 

While the PCAC masses renormalize multiplicatively, the bare Lagrangian 
masses require an additional additive renormalization. In the difference d=m^ are - 
mjj™, this additive renormalization is eliminated. Moreover, the multiplica- 
tive renormalization factors cancel in the ratio r=m^ CAC /m^ AC . To obtain 
fully renormalized quantities, we must still multiply d by 1/Z$, the inverse 
of the scalar density renormalization factor. From the renormalized mass dif- 
ference d/Zs and the renormalization independent ratio r we obtain m™£ — 
(d/Zs)/(r — 1) and m™ 11 = (rd/Zs)/(r — 1) in the unimproved case. Our final 
analysis is tree- level 0{a) improved with slightly more complicated formulae 
(see Sec. 11.2 of (H). 

To compute Zs nonperturbatively in the RI scheme, we apply the Rome- 
Southampton (RS) method 3J| with tree-level improvement, augmented with 
nonperturbative running. Our procedure eliminates the possible difficulties of 
the RS method on coarser lattices. Since the RI scheme is defined in the Nf—3 
(i.e. with three degenerate quarks) chiral limit, we generate Nf=3 configura- 
tions at our five lattice spacings and, for each /3, at four or more values of m q 
to allow an extrapolation to the massless limit. Thus for each /3 we compute 
the renormalization constant Z^ l (f3, /i), at renormalization scale /i, as described 
in [32|]. The RS procedure defines a valid renormalization scheme as long as 
/i<^.ir/a. However, only if /i3>Aqcd can the results be converted perturbatively 
to other schemes (including intrinsically perturbative schemes such as MS) or be 
used in a perturbative context. On coarser lattices, it is difficult to simultane- 
ously satisfy both constraints on ji. To solve this difficulty we first determine the 
quark masses at /1 = 1.3 and 2.1 GeV, then apply continuum nonperturbative 
running, as defined in [35j . up to [i! = 4 GeV. 

(iv) Combined analysis of mass and lattice spacing dependence. For the 
masses, two strategies, called "Taylor fit" and "chiral fit" [l[ are applied. Clearly, 
the results of these fits are dominated by the results at the physical point. In the 
analysis, two different pion mass ranges are used, namely M T <340, 380 MeV. 

The strange and average up-down quark masses renormalized in the RI 
scheme at 4 GeV are extrapolated to the continuum and interpolated to the 
physical mass point. In these fits, we include terms to correct linear (a s a) or 
quadratic (a 2 ) effects. A combined mass and lattice spacing fit is carried out. 
We show the continuum extrapolation for m uc i and m s in the RI scheme at 
4 GeV, as well as their ratio, in Figure [5J In order to control the systematic 
uncertainties we carry out 288 such analyses 32]. The figure depicts results 
from one analysis with one of the best fit qualities. 

Our procedure yields the RI quark masses m u d and m s , with statistical and 
fully controlled systematic errors. These results do not rely on perturbation 
theory and from them it is straightforward to obtain the quark masses in other 
commonly used frameworks such as renormalization group invariant (RGI) and 
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Figure 2: Continuum extrapolation of the average up/down quark mass, of the strange quark 
mass and of the ratio of the two. The errors of the individual points, which are statistical only 
here, are smaller than the symbols in most of the cases. The only exceptions are the light 
quark mass and its ratio to the strange quark mass at the two finest lattice spacings. These 
exceptions underline the importance of using physical quark masses to reach a high accuracy. 





RI(4GeV) 


RGI 


MS(2GcV) 


m s 
m ud 


96.4(1. 1)(1.5) 
3.503(48)(49) 


127.3(1.5)(1.9) 
4.624(63)(64) 


95.5(1.1)(1.5) 
3.469(47) (48) 


m u 
md 


2.17(04)(10) 
4.84(07)(12) 


2.86(05)(13) 
6.39(09)(15) 


2.15(03)(10) 
4.79(07)(12) 



Table 1: Renormalized quark masses in the RI scheme at 4 GeV, and after conversion to RGI 
and the MS scheme at 2 GeV. The RI values are fully nonperturbative, so the first column is 
our main result. The first two rows emerge directly from our lattice calculation. The last two 
include additional dispersive information. 



MS [36( ones. Moreover in [32j, we show that 4 GeV in the RI scheme is suffi- 
ciently large that the perturbative running required to obtain RGI masses and 
the matching to the MS scheme at 2 GeV both yield subdominant uncertainties. 

The determination of the individual up and down quark masses at the phys- 
ical point is in principle possible using exclusively lattice simulations. To that 
end one may include the electromagnetic U(l) gauge field into the lattice frame- 
work, as was done recently in [37[. Such a project goes beyond the scope of the 
present paper, which deals with QCD only. Nevertheless our precise m s and m u d 
values can be combined [32| with model-independent results based on dispersive 
studies of rj^Sw decays to derive the individual up and down quark masses (c.f. 
Tab. [I}. In this approach the relationship between the input parameters and 
experiments is not as transparent as for the determination of m s and m u d (see 
condition 6 above). 

We have performed a lattice QCD determination of the light quark masses 
which includes all of the ingredients that we believe are required to achieve full 
control over systematic errors. In particular, we have eliminated the need for 
difficult chiral extrapolations by performing simulations all the way down to the 
physical pion mass (and even below); gained control over the necessary contin- 
uum extrapolation by performing simulations at five lattice spacings down to 
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a = 0.054 fm; eliminated perturbative uncertainties by performing the renormal- 
ization and running fully nonperturbatively; and controlled the infinite volume 
extrapolation by working with lattice sizes up to 6 fm. 

The precision reached for m s and m uc i is somewhat below the 2% level and for 
m s /m„d=27.53(20)(08), which is scheme-independent, it is better than 1%. For 
m u and m^, which include additional dispersive information, it is about 5% and 
3%, respectively. Despite their use of significantly different methods, MILC (l7j . 



RBC [21[ and HPQCD [191 ] . the three collaborations which have performed the 
most extensive 7V/=2+l computations besides ours, obtain results for these 
masses which are at most 1.5 combined standard deviations away from ours. 

Our results provide precise and reliable input for phenomcnological calcula- 
tions which require light quark mass values. They highlight the progress that 
has been achieved since the early days of quark mass determinations [38| by 
showing that phenomenologically relevant lattice QCD calculations can now be 
carried out bracketing the physical values of the light quark masses. 

The details of this work can be found in _32j . 
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